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ABSTRACT 



Context. The knowledge of mass loss rates due to thermal winds in cool dwarfs is of crucial importance for modeling evolution of 
physical parameters of main sequence single and binary stars. Very few, sometimes contradictory, measurements of such mass loss 
rates exist up to now. 

Aims. We present a new, independent method of measuring an amount of mass lost by a star during its past life. 
Methods. It is based on the comparison of the present mass distribution of solar type stars in an open cluster with the calculated 
distribution under an assumption that stars with masses lower than M lim have lost an amount of mass equal to AM. The actual value 
of AM or its upper limit is found from the best fit. 

Results. Analysis of four clusters: Pleiades, NGC 6996, Hyades and Praesepe gave upper limits for AM in three of them and the 
inconclusive result for Pleiades. The most restrictive limit was obtained for Praesepe indicating that the average mass loss rate of cool 
dwarfs in this cluster was lower than 6 x 10 _11 Mo/yr. With more accurate mass determinations of the solar type members of selected 
open clusters, including those of spectral type K, the method will provide more stringent limits for mass loss of cool dwarfs. 

Key words, stars: activity - stars: late-type - stars: mass loss - stars: winds, outflows 
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1. Introduction 

Chromospheric-coronal activity is ubiquitous in cool stars pos- 
sessing subphotospheric convection zones. Hot coronas result in 
existence of thermal winds carrying away stellar mass and an- 
gular momentum. The precise knowledge of the mass loss rate 
(MLR) during the stellar life is of crucial importance for mod- 
eling evolution of mass, luminosity and rotation of single stars 
and of orbital parameters of close binary stars. Mass loss of a 
central star of a planetary system may also influence evolution 
of physical parameters of orbiting planets. Proximity of the Sun 
made possible direct measurements of solar wind with the use of 
interplanetary probes. The average velocity of the wind near the 
Earth orbit is of the order of 400 km s _1 and the average MLR 
is 2 x lO~ 14 M0/yr with uncertainty of some 50 % (Feldman et 
al. 1977). The expected MLRs for other cool dwarfs may dif- 
fer from that at most by a couple of orders of magnitude but 
even the highest MLRs of such stars are far lower than those 
occurring in hot, early type stars or red giants. The present ob- 
servational capabilities are still insufficient to obtain the data 
for other dwarfs with comparable to the solar case accuracy. 
For many years only indirect methods of MLR estimate of such 
stars existed. Analyzing the infrared and radio observations of 
M dwarfs Mullan et al. (1992) obtained a value of lO" lo M0/yr 
as an upper limit for MLR of these stars. Later, van den Oord & 
Doyle (1997) revised that limit down to lO~ 12 M0/yr. A similar 
value was obtained by Lim & White (1996). Wargelin & Drake 
(2002) obtained an upper limit of 3 x 10~ I3 MQ/yr for Proximo 
Cen. In the recent years Wood and coworkers published a num- 
ber of papers in which they analyzed profiles of L a arising at 
the collisional front between stellar wind and interstellar mat- 



ter. By modeling the collision process the authors were able to 
measure MLRs for a dozen stars (Wood et al. 2002, 2005). The 
measured rates extend from 2 x lO~ 12 M0/yr for 70 Oph down 
to 3 x lO~ 15 M0/yr for DK UMa. The authors derived a relation 
between MLR and age (or X-ray flux). Their latest results give 
M oc r 233 for t > 0.7 Gyr but observations of younger, more 
active stars suggest a sudden drop of MLR from about 2 x 10~ 12 
to ~ 10- 13 M o /yr, i. e. by more than an order of magnitude. 
This unexpected result is based on scanty data and needs con- 
firmation. A different, indirect method of estimation an amount 
of mass lost by the Sun during its past life was suggested by 
Sackmann & Boothroyd (2003). Based on the geological evi- 
dence that the terrestrial atmosphere was warm in early life of 
the Solar System and liquid water was present on Mars, they cal- 
culated the minimum initial mass of the Sun needed to produce 
enough luminosity, and reconstructed time variation of the so- 
lar MLR so adjusted that the favorable thermal conditions on the 
Earth and Mars have been kept as the Sun aged. The optimum 
conditions were reproduced for the initial mass of 1 .07 Mq and 
MLR decreasing exponentially with time. 
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Here we present another method of measuring the total 
amount of mass lost during the stellar life. It is based on the 
analysis of the stellar mass distribution in an open cluster. The 
next section describes details of the suggested method. The last 
section contains discussion of the results and conclusions. 
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2. Method 



2.2. Method of the mass loss estimate 



2.1. General considerations 

The total energy flux carried away with the stellar wind, F t0( , 
consists of three components (Holzer 1987) 



Ffot — Foray ~t~ Fqo ~t" Fn 



(1) 



where F srdv is the energy flux needed to carry the wind matter 
out of the potential well of the star, F «, is the kinetic energy 
flux of the wind in infinity and F ra( j is the flux radiated away by 
the wind. Based on observations of the solar wind F m & can be 
neglected. We thus obtain 



f tot = \M(vl c + vt) 



Mv 2 esc . 



(2) 



where v eS c is the escape velocity from the stellar surface and 
Voo is the wind velocity in infinity. The last approximate equality 
results from the assumption that both velocities are of the same 
order. With v eS c = (2GM/R) 1/2 , where G is gravitational con- 
stant, and M and R are stellar mass and radius, respectively, we 
have 



2GMM 
R 



(3) 



The ultimate source of energy for thermal winds in cool 
stars is stellar luminosity which drives convection and, via a 
chain of physical processes, makes a magnetized wind to blow. 
Rotational energy, albeit necessary for magnetic field generation 
(dynamo does not work in nonrotating stars), is of secondary im- 
portance. If so, the total energy flux of the wind can be expressed 
as a fraction kL of stellar luminosity 



M - 1.5 x 10 



-8 



kLR 
~M 



(4) 



where M, R and L are in solar units and M is expressed in 
solar mass per year. The observed solar value of MLR gives 
&0 = 1.3 x 10~ 6 . For comparison, Kudritzki & Reimers (1978) 
derived an empirical relation for M-type giants 



i -i LR 

M = 5.5xl(T 13 — , 
M 



(5) 



from which £mg = 3.7 x 10~ 5 . Surprisingly, although MLRs 
of M-giants are some 7 orders of magnitude higher than the solar 
value, the efficiency factor &mg, describing the transformation of 
the red giant radiation into wind energy, is only about 30 times 
higher than kQ. 

How high values can MLR reach in the most active cool 
dwarfs? Stars with X-ray emission at the saturated level have 
L, « 10~ 3 L and if we apply the same efficiency factor to MLR 
(i. e. k — 10~ 3 ) we obtain an upper limit of 1 .5 x 10~' 'MQ/yr for 
a one solar mass star. The upper limit drops, however, sharply 
for lower mass stars due to a fast decrease of their luminosity 
and reaches a value of 7 x 10~ 13 MQ/yr for a 0.4 solar mass star. 
(Note that, according to the results of Wood et al. 2005, the pro- 
portionality between MLR and X-ray luminosity breaks down 
for the most active stars). This brief discussion indicates that the 
expected MLR of active solar type stars should not considerably 
exceed ~ 10 _1 'Mo/yr and yet it may be significantly lower if the 
results of Wood et al. (2005) are confirmed. 



Initial mass distribution of a stellar cluster is well described by 
the Salpeter function (Salpeter 1955). Recent discussion of this 
distribution by Kroupa (2002) showed, as the author stresses, 
that the initial mass function of different populations of stars 
shows an extraordinary uniformity and is well described by the 
Salpeter function, contrary to the simple minded reasoning that 
it should vary with star-forming conditions. The Salpeter mass 
distribution is given by 



AM 



(6) 



where AN is a number of stars with masses between M and 
M + AM. Values of a determined for individual star aggregates 
in the intermediate mass range concentrate around the original 
Salpeter value a = 2.35 with a scatter expected from statistics 
(Elmegreen 2001). 

The mass distribution evolves with the cluster age due to 
different dynamical effects, like evaporation of low-mass stars 
and stellar evolution affecting mostly high-mass stars. The mass 
distribution of solar type stars (with masses, say, between 0.5 
and 2 Mq) is expected to be least affected by the above men- 
tioned effects. It will, however, evolve due to mass loss via ther- 
mal wind from stars with subphotospheric convection zones. Let 
us assume that the mass distribution of solar type stars with 
M > M\{ m has not changed during the past cluster life but stars 
with M < My lm have lost an amount of mass equal to AM. The 
assumption is justified by the sharp rise of coronal activity at 
mid-F spectral type and, equally sharp, drop of stellar rotation 
velocity around the same type, indicating a sudden appearance of 
magnetized winds over a narrow spectral range (Schmitt 1997; 
Gray 2005). Let AM be independent of mass. A gap in the mass 
distribution then occurs just below My lm with a width of AM and, 
at the same time, the distribution function for lower mass stars 
will run lower than the initial function. Fig Q] shows the initial 
(Salpeter) mass distribution and two distributions modified by 
mass loss of 0.1 and 0.3 Mq when My lm = 1.25Mq. The second 
value of AM is probably unrealistically high but it is shown here 
to emphasize the expected distribution modifications. 

There are indications that the bulk of mass loss occurs when 
stars are young and active i. e. younger than about 1 Gyr. The 
time-scale of spin-down of young, rapidly rotating stars is about 
1 - 2 x 10 8 years (Stepien 1988) which means that after several 
hundred Myr the activity drops considerably. The observations 
of stellar X-ray activity by Giidel, Guinan & Skinner (1997) 
show, indeed, that the activity decreases by about 1.5 order of 
magnitude during the first Gyr of stellar life. These results indi- 
cate that the optimum cluster age for the analysis of AM is from 
several hundred Myr, up to about 1 Gyr. Stars in very young 
clusters have not had enough time to lose an appreciable amount 
of mass whereas AM in old clusters is expected to increment 
negligibly with age due to a decrease of MLR down to the solar 
value or even beyond. At the same time, mass distribution of an 
old cluster becomes notably perturbed by other effects, like evap- 
oration of low mass stars and evolution of massive stars. With 
AM known, the average MLR can be estimated: M av = AM/T, 
where T is the cluster age. 

Unfortunately, accurate mass values of solar mass mem- 
bers are presently known only for very few clusters in the op- 
timum age range. We selected four clusters for the analysis: 
Hyades, Praesepe, Pleiades and NGC 6996. Masses of Hiads 
were determined by Perryman et al. (1998). Later analysis by 
de Bruijne, Hoogerwerf & de Zeeuw (2001) confirmed these 
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Fig. 1. Expected differences among stellar mass distributions in a 
cluster, resulting from the assumed mass loss: solid line, labelled 
dO - no mass loss; dotted line, labelled d-0.1 - mass loss of 0.1 
Mq in stars with M < 1.25M©; broken line, labelled d-0.3 - 
mass loss of 0.3 Mq in the same stars. All three distributions 
have a common part for stars with M > 1 25Mq (the vertical 
dotted and dashed lines at M = 1 25Mq are slightly displaced 
from their exact positions for better visibility). Note a gap in two 
latter distributions, broadening with the increased mass loss. 

values. Perryman et al. (1998) give masses for 218 stars. After 
rejecting all giants, variables and stars marked by de Bruijne et 
al. (2001) as nonmembers, 133 dwarfs were left. Data on stel- 
lar masses of Praesepe and Pleiades were taken from Raboud 
& Mermilliod (1998a, b). The authors selected 185 members 
of Praesepe and 270 members of Pleiades and determined their 
masses from colors using isochrones with metallicity Z = 0.02. 
Masses of known binary components were determined as de- 
scribed by Raboud & Mermilliod (1998a). No published masses 
for members of NGC 6996 exist so we computed them using 
UBV photometry of Villanova et al. (2004). The following val- 
ues of the parameters, needed to convert the observations into 
absolute photometry, were adopted: distance moduli 9.4 mag, 
E(B - V) = 0.54 and A v = 3E(B - V). Bolometric correc- 
tions were applied following Bessel, Castelli & Plez (1998). 
Individual masses were calculated from the relation (Lang 1992) 

1Og (^) = - 2631Og fe)- a021 - (?) 

Photometric errors result in +0.05Mq random mass errors. 
The errors of absolute values of stellar masses are likely to be 
larger, due to probable systematic errors, but the uncertainties 
of the relative mass determinations should be close to those re- 
sulting from photometry. Errors of mass determination of sin- 
gle stars in the other analyzed clusters are assumed to be of 
the same order i. e. +0.05Mq. In case of resolved binaries with 
no spectroscopic orbits expected errors are several times larger. 
Unresolved binaries present another problem. They enter the 
analysis as single objects instead of two stars, producing thus 
a bias in the observed mass distribution. Fortunately, the com- 
monly employed methods of detecting binaries are most sensi- 
tive to systems with similar components, i.e. when mass ratio q 
is close to 1 . As a result, frequency of known binaries among so- 
lar type stars in the investigated clusters is as high as, or higher 
than, among nearby field stars for which binary searches are be- 



Table 1. Results of a comparison of the calculated and observed 
mass distributions. The compared distributions are identical at 
the significance level of 5% if statistic A < A\i m = 1.36. 



Cluster 


AM 


A 


Cluster 


AM 


A 


Pleiades 





2.18 


Hyades 





0.58 




0.1 


2.96 




0.1 


0.93 




0.2 


3.81 




0.2 


1.42 


NGC 6996 





0.55 


Praesepe 





1.08 




0.1 


1.05 




0.1 


1.75 




0.2 


1.58 




0.2 


2.62 



lieved to be nearly complete (Halbwachs et al. 2004). This in- 
dicates that in well studied clusters the percentage of missing 
binaries with both component masses from the considered mass 
range [0.5,2.0] Mq is low, except, possibly, in the cluster cen- 
tral regions where crowding of stellar images makes analysis 
more difficult. A stronger bias is expected among low q bina- 
ries with a solar type primaries and very low mass secondaries. 
Only primaries enter the statistics and secondaries are missing, 
which produces a deficit in the lowest mass tail of mass distribu- 
tion (Kroupa 1995). Such a deficit does not matter in our case as 
long as we restrict the analysis to solar type stars only. 

We assume conservatively that two mass distributions differ- 
ing in AM by less than 0.05Mq are statistically indistinguish- 
able, so if for example the original Salpeter function correctly 
describes the observed distribution of a given cluster, we adopt 
0.05 Mq as an upper limit for mass loss of its solar type mem- 
bers. 

The observed mass distribution in each analyzed cluster was 
binned in 0.1 Mq wide intervals. Multiple star components 
with known masses were treated as single stars. The calculated 
distribution, corresponding to a given value of mass loss, was 
calibrated against the mass interval 1.25 - 2.0Mq and its part 
corresponding to M < \.25Mq was compared to the observa- 
tions. The Kolmogorov-Smirnov test was applied to verify a null 
hypothesis that the calculated distribution is identical with the 
observed one in this mass interval. The test uses the statistic 
A (sometimes denoted by D„) describing the maximum differ- 
ence between the empirical and the given (theoretical) cumula- 
tive distribution functions. We adopted the significance level of 
5% which corresponds to the statistic A < Aum = 1-36. The pro- 
cedure was applied to the Salpeter function and the distributions 
with a few values of AM. The results are shown in Table [T] 

3. Results 

3.1. Comparison of the calculated mass distributions to 
observations 

Assuming that solar type members of a given cluster have lost 
AM C \ mass over their past life, we expect the function A(AM) 
to reach a minimum A - /l m ; n at AM = AM C \. If the condi- 
tion /Imin < Aiim is fulfilled we assume that AM C \ is the most 
probable value of the amount of mass lost by the cluster mem- 
bers although any value of AM from the neighbourhood of AM C \ 
is also possible as long as A(AM) < A\\ m . The latter condition 
determines the uncertainty of the obtained mass loss. Its rigid 
application may, however, result in an underestimation of actual 
uncertainty of AM C \ when A^ is only barely lower than Ay lm and 
A(AM) quickly exceeds A\i m for AM moving away from AM C \. To 
obtain a more realistic measure of the actual accuracy of AM Q \ 
both, a depth and a shape of the A minimum should be consid- 
ered. A deep, sharp minimum with A rising rapidly away from it 
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means that AM C \ is well determined with a high degree of con- 
fidence. A broad, shallow A minimum means that the value of 
AM c i is poorly constrained. 

TableQ]shows that /t m ; n is reached for AM = for all the in- 
vestigated clusters, which means that only an upper limit for the 
mass loss of each cluster can be determined. In case of Pleiades 
- the youngest of the discussed clusters (see Table |2j, even the 
distribution with AM = poorly describes the observed distri- 
bution. This may be connected with an apparent deficit of mas- 
sive stars among the presently observed members (Moraux et al. 
2004). Fig. [2] shows the observed distributions of all four clus- 
ters with the Salpeter function overplotted. For NGC 6996 and 
Hyades the distributions with AM = 0.1 Mq also describe satis- 
factorily the observed distributions at the assumed significance 
level so the upper limit for mass loss in these clusters is equal 
to 0.15 Mq. Note, however, that in case of Hyades the rise of 
A(AM) is so slow that a value of 0.20 Mq only weakly violates 
the 5 % condition, so the limit of 0.25 Mq is nearly equally 
probable as 0.15 Mq. The same limit of 0.25 Mq was obtained 
from a comparison of the calculated distributions to observations 
of 90 single Hyads with best known masses, carefully selected 
by de Bruijne et al. (2001). Finally, only the Salpeter function 
describes satisfactorily the mass distribution in Praesepe and the 
rise of A(AM) is quite sharp (Table [TJ. While the probability of 
occuring A > 1.08 is about 0.20, it drops to less than 0.01 for 
A > 1.75. We adopt a value of 0.05 Mo as an upper limit for 
mass loss in this cluster. Dividing the upper limit by the clus- 
ter age we obtain an upper limit for an average MLR over the 
past cluster life. The results (together with the adopted ages) are 
shown in Table [2] 

As we see, the derived MLRs in NGC 6996 and Hyades 
are not very restrictive. However, the upper limit for mass loss 
and MLR in Praesepe look interesting. The maximum mass loss 
is lower than the solar mass loss of 0.07 Mq determined by 
Sackmann & Boothroyd (2003). The maximum MLR is close to 
the theoretically estimated upper limit (see above) and is lower 
than some of the empirical determinations (e. g. Mullan et al. 
1992). 

The result for Praesepe demonstrates a potential usefulness 
of the suggested method. The presently attained accuracy of 
the method is limited by two main factors: the accuracy of 
mass determination and scarcity of data, particularly for lower 
mass members. The present mass determinations extend only to 
M ~ O.8M0. Deeper surveys, including stars of spectral type 
K, i. e. down to 0.5 MQ, are needed, together with an increased 
accuracy of mass determination. When masses of stars from the 
spectral range A-K in stellar clusters with age ~ 1 Gyr are known 
with relative accuracy of 0.01-0.02 Mq, it will be possible to de- 
tect average MLRs at the level of ~ 10 _1 1 Mq or less. Such an 
accuracy should be possible to reach with new generation in- 
struments, like LSST, Gaia and fibre-fed multiobjects spectro- 
graphs. The residual unresolved binaries will add to the obser- 
vational noise smearing the notch in the predicted mass distribu- 
tion (Fig.Q]). Fortunately, the Kolmogorov-Smirnov test is only 
weakly sensitive to such an effect. 

Systematic errors may be larger but this does not hurt - such 
errors influence only absolute mass determinations which are ir- 
relevant here. Moreover, they can be offset by adequately vary- 
ing the value of the limiting mass for mass loss to occur. If the 
future data are accurate enough to detect a possible change of the 
slope of the mass distribution with age it may be possible to an- 
alyze mass dependence of MLR. This, however, needs a further 
increase of accuracy of mass determination and an extension of 
the observed mass distribution to still lower masses. 




1.0 1.2 1.4 1.6 1.8 3.0 




Fig. 2. The observed mass distributions in three clusters with the 
Salpeter functions overplotted. The Salpeter function is normal- 
ized in each case to the distribution of stars with masses from 
the interval 1.25-2.0 Mq. 

Table 2. Upper limits for mass loss and the resulting upper limits 
for average MLRs of stars in the investigated clusters. 



Cluster 


Age 


Max. mass loss 


Max. average MLR 




Myr 


Mq 


10- u M o /yr 


Pleiades 


100 


? 


? 


NGC 6996 


325 


0.15 


43 


Hyades 


625 


0.15-0.25 


24-40 


Praesepe 


832 


0.05 


6 



3.2. Conclusions 

A new method of mass loss determination in cool dwarfs is pre- 
sented. It is based on a comparison of the observed mass dis- 
tribution in a stellar cluster with a distribution modified by an 
assumed mass loss AM occurring during the cluster past life. 
The method was applied to four open clusters with ages between 
100 Myr and 1 Gyr. The most significant result was obtained 
for Praesepe for which an upper limit of 0.05 Mq for AM was 
determined. From the age of the cluster an upper limit for aver- 
age MLR can be calculated. It is equal to 6 x 10~ n Mo/yr. This 
is a promising result. With more complete surveys of cluster 
members and better mass determinations, a significantly better 
estimate will be obtained. 
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